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We discuss some aspects related to the so-called Hilbert space Average Method, as an alternative 
to describe the dynamics of open quantum systems. First we present a derivation of the method 
which does not make use of the algebra satisfied by the operators involved in the dynamics, and 
extend the method to systems subject to a Hamiltonian that changes with time. Next we examine the 
performance of the adiabatic quantum search algorithm with a particular model for the environment. 
We relate our results to the criteria discussed in the literature for the validity of the above-mentioned 
method for similar environments. 



INTRODUCTION 



The study of open quantum systems has attracted re- 
newed attention during the last years. One important 
reason for this is the expected advent of future quantum 
computers [l| . The interaction of the quantum computer 
with its surroundings can introduce some degree of de- 
coherence which can, eventually, ruin the performance of 
the quantum algorithm. Although this phenomenon can 
be mitigated with the help of error-correction methods, 
a deeper understanding of how the ambiance operates on 
the smaller system can also be used to improve the work- 
ing conditions. In fact, recent papers have shown that, 
for some models of system-ambiance interaction, the loss 
of coherence can be smaller if the ambiance temperature 
is increased @, 0|- Similarly, designing some engineered 
reservoirs with controlled coupling and state of the envi- 
ronment can reduce the decoherence rate @, One can 
even consider purely dissipative processes, which turn out 
to be equivalent to a quantum circuit model for quantum 
computation 0|. On the other hand, systems subject to 
decoherence will experience a transition from a quantum 
to a classical state. The study of this transition will give 
more insight about the nature of Quantum Mechanics 
and its differences with a classical perception [3]. 

There are different approaches which have been devel- 
oped in the literature in order to describe the evolution 
of open systems, based on different techniques such as 
master equations or superoperators [§]. As an alterna- 
tive to these methods, we will study the behavior of the 
open system using the so-called Hilbert space Average 
Method (HAM, in what follows) @, [ijj]. This method 
has been proved to give, in some situations, better results 
than conventional Time Convolutionless (TCL) approxi- 
mations, and comparable to correlated projection super- 
operator techniques . We will extend this approach to 
the case of a time- varying Hamiltonian acting on the open 
system, and will show that, under a suitable choice of the 
operators defining the HAM scheme, the resulting equa- 
tions will adopt the form of traditional master equations, 
at least up to second order in the system-environment 
coupling. 



As an application, we will consider the case of a 
quantum system which is designed to perform a Grover 
search FlJ EE E3| via adiabatic quantum computation 
nam, 1 1 71. Our purpose is to analyze the response 
of the quantum computer when coupled to the external 
influence of an environment, introduced with the help of 
some specific model. This problem has been considered 
by several authors @, [H, M, HI H3, 0, HI EE H3| , 
but here we relate it to similar models for the bath 
which have discussed within the HAM formalism. We 
compare our results with a numerical simulation of the 
Schrodinger equation obeyed by the full system. 

In Section II we introduce some basic notations. In 
Sect. Ill we briefly revisit the HAM method, using an 
approach that makes not use of the algebra of the in- 
volved operators. An approximated evolution equation, 
extended to the case of a time-changing Hamiltonian act- 
ing on the system, is obtained in Sect. IV, and we also 
make a connection with familiar master equations. Sect. 
V is devoted to the analysis of adiabatic search when the 
quantum computer interacts with a particular environ- 
ment. The evolution of the system is followed both by 
the exact Schrodinger equation and by solving the ob- 
tained approximated equations. The comparison of both 
calculations is discussed within the framework of known 
criteria for similar models, which have derived within the 
HAM formalism. Our results are summarized on Sect. 
VI. 

We work in units such that h = 1 . 



II. BASIC NOTATIONS 

We wish to study the evolution of an open system (S) 
in contact with an environment (E). First we introduce 
the basic quantities in the Schrodinger picture, and then 
we will define a convenient interaction picture for this 
problem. System S is subject to a time-dependent Hamil- 
tonian Hs(t) . Let us denote by He the Hamiltonian de- 
scribing the free evolution of the environment, and by V 
the interaction between both systems. Both He and V 
are assumed to be time-independent. Moreover, we make 
the hypothesis that [Hs(t), He] = . The evolution of 
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the density matrix p{t) of the complete (S+E) system is 
therefore given by 

j t p{t) = -i[H{t),p{t)] (1) 

with H{t) = H s (t) + H E + V the total Hamiltonian. 
We can define an interaction picture density matrix pi(t) 
as follows: 

Pl {t) = exp(iH B t)p(t) expHfffii). (2) 
The equation for pi(t) is easily obtained: 

j fPl (t) = -i[H s (t), Pl (t)} - i[VK*), />/(*)]. (3) 
where 

Vi(t) = exp(iH E t)V exp(-iH E t) (4) 

is the interaction operator in the interaction picture. In 
what follows we will assume, unless otherwise specified, 
that we are working in the above-defined picture, and 
will therefore omit the subscript T. 



III. HILBERT SPACE AVERAGE METHOD 

In this section we briefly review the Hilbert space Av- 
erage Method (HAM), as an alternative to describe the 
dynamics of an open quantum system. A more detailed 
description of the method can be found in 0, [n], EH ■ 

The idea is to replace the dynamics of the density ma- 
trix p(t) describing the full system (i.e. open system plus 
reservoir) by an effective density matrix a{t) which is 
simpler to describe, with the condition that the expected 
values of a given set of operators is reproduced. Let us 
assume that we are interested on a set of operators {P n } 
and we want to define a density matrix satisfying the 
boundary conditions 

Tr[a(t)P n ] = p n (t), (5) 

where the functions p n (t) are assumed to be known 
(actually, they will be determined by the dynamics), and 
Tr stands for the trace over the whole Hilbert space. We 
would like to determine a(t) respecting the above con- 
ditions, but otherwise unknown. To this end we estab- 
lish the following procedure. We maximize the entropy, 
in order to account for our ignorance about the effec- 
tive density matrix, but add constraints corresponding 
to Eqs. under the form of Lagrange multipliers. For 
our purpose, it is simpler to consider the linear entropy 
S[p] = 1 — Tr[p 2 ]. In this way, we find the extrema of 
the functional 

I[a] = S[a] - a nTr[aP n ] (6) 

n 



with {a n } the above defined Lagrange multipliers. Fol- 
lowing this procedure, one arrives to the expression 

a(t) = J2bn(t)P n , (7) 

n 

with b n — —a n /2 , and we have made explicit the 
time dependence of a(t) . The new functions b n (t) are 
determined by Eqs Of course, one also has to make 
sure that the condition Tr[a(t)] = 1 is satisfied. 

In [{| [l(| EH, the authors introduce the HAM method 
as an average over all possible states in the Hilbert space 
that accounts for conditions Eqs. © (here is where the 
name HAM comes from) , and introduce the set of oper- 
ators {P n } as obeying a closed algebra. In contrast, our 
derivation of Eq. ([7]), although perhaps less intuitive, 
makes no use of the algebra of the operators {P n } ■ 



IV. APPROXIMATED EVOLUTION EQUATION 

The exact dynamics of a(t) will be given by solving 
an equation like Eq. |(3]) in the interaction picture. This 
equation does not admit a simple, closed form, solution. 
In this section we will investigate some approximation 
that is easier to solve, and will allow us to make a connec- 
tion with standard master equations. We first introduce 
the evolution operator U(t + r, t) from instant t to t + r. 
We then have 

a(t + t) = U(t + t, t)a{t)U\t + r, t). (8) 

In order to separate the evolution due to H s {t) from 
that due to V{t) we consider a sufficiently small r and 
approximate U(t + r,t) by a second-order Suzuki decom- 
position [H] 

U(t + r, t) ~ exp(-^H s (t))D v (t + T,t) exp{-i^H s {t)), 

(9) 

where Dy(t + T,t) is the evolution operator describ- 
ing the time evolution due to V(t) alone (i.e., neglecting 
Hs(t) in the Hamiltonian), and verifies the equation 

i-D v (t,t ) = V{t)D v (t,t ) (10) 

Inserting Eq. |(9]) into (JSj) and expanding the exponentials 
in powers of r gives 

a{t + r) = a(t) + Aa(t + r,t) + ir[a(t), H s (t)] + C(r 2 ), 

(11) 

with the definition 

Aa(t + r,t) = D v (t + T,t)a(t)D\r(t + T,t) - a(t). (12) 

We now would like to make an approximate treatment 
of the quantity Aa(t + r, t) defined in the latter equation. 
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To this purpose, we use the solution of Eq. (fTO)) up to 
second order in the potential 

/t+T pt+T pS 

dsV(s)- / ds ds'V(s)V(s'). 

(13) 

Within this approximation, Eq. Il"12|) reads 

[■t + T 

Aa(t + r,t) = -i J ds[V(s),a(t)} 

- ds j* ds'[V(s), [Vis'), a(t)}] . (14) 

Starting from this approximation, one can derive the cor- 
responding equations for the functions {p n {t)} , following 
the procedure described in [lflj. To this end, one needs 
to specify the operators {P n } and the algebra verified by 
them. 

We can also establish a connection with familiar master 
equations, which is done in a trivial way within the above 
formalism. We simply assume that the effective density 
matrix a(t) can be factorized as 

a(t) = p s (t) <g> p E (15) 

where is a density matrix that approximates the 
state of the bath, and ps{t) is the density matrix for the 
system S, related to a(t) via 

p s (t)=Tr E [a(t)], (16) 

and Tte indicates the partial trace over the environ- 
ment E. 

Let us introduce an orthonormal basis {\i >} in the 
Hilbert space corresponding to system S, and write Eq. 
(fT5l) in the following way: 

a{t) = Y J P 3l { t )^><j\®PE (17) 

By comparing with Eq. J7]) we identify the operators 
P n associated with the ansatz ffl5|) 

Pn = \i><j\®PE (18) 

where n indicates a given pair i,j. One also easily 
obtains from Eq. ((5j) that 

p n {t) = Pii{t)Tr[p%]. (19) 

Notice that the functions Pij (t) are related to the ma- 
trix elements psi,j(t) in the basis {\i >} via Pij{t) — 
PSji(t) ■ 

In order to obtain a more detailed expression for the 
quantity Aa(t + r,t) one needs to specify the interaction 
V . Let us assume that this operator is defined, in the 
Schrodinger picture, by 

V = ^Ai®d (20) 



with the Hermitian operators Ai (C,) acting on the 
system S (E). Of course, one can consider a more general 
situation where these operators are not Hermitian, and 
simply add the Hermitian conjugate to V . However, a 
simplified version like Eq. (|20|) will be sufficient to our 
purposes. In the interaction picture, the above formula 
becomes 

V I {t)=Y,A i ®C i (t) (21) 

i 

with Ci{t) = e?q>(iH st)Ci e~>q>(—iH Et) . Hereafter, we 
omit the subindex T, as anticipated in Sect. II, and 
assume that we are working in the interaction picture. 

In what follows, we are interested in the difference 

Ap s (t+T, t) = ps{t+r)-p s (t) = Tr E [Aa(t+T, t)}. (22) 

We will obtain an approximation to this quantity by 
using Eqs. (fT4|l and l|2ip . The rest of this section is a 
standard manipulation which is common to the deriva- 
tion of master equations. Our purpose is only to show 
that the program we started in Sect. Ill does indeed lead 
to such kind of equations. The interested reader is ad- 
dressed to the existing bibliography (see, e.g. 3). The 
final expression reads 

j t p s (t) = ~i[H s (t),p s (t)} 

+ ?ik{t){A k p s {t)Ai - AiA kPs (t)} + h.c. (23) 

l,k 

where the Hermitian conjugate refers only to the sum- 
mation. In obtaining Eq. (|23|l we have taken the limit 
t — > , and we have defined 

T lk {t) = lim - [ ds f ds'G lk {s,s'), (24) 

with 

G ik {s, s') = Tr E {Ci{s)C k {s')p E } (25) 

the bath correlation functions. We also have made the 
usual hypothesis that 

Tr E {C l {s)p E } = 0. (26) 

Eq. f23|) takes then the familiar form of a master equa- 
tion, which becomes of the Lindblad type in the case that 
the coefficients Ti k (t) are independent of time. 

V. MODEL FOR ADIABATIC SEARCH AND 
INTERACTION WITH THE ENVIRONMENT. 

We now analyze a particular and interesting exam- 
ple, which can be cast either under the form of HAM 
or master equations, according to the discussion of the 
previous section. We study the performance of an 
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open quantum system (the quantum computer) consist- 
ing on n qubits, while it does an adiabatic search for 
a marked state |m > out of N = 2" possible con- 
figurations, subject to the interaction with an environ- 
ment. This problem has been addressed by several au- 
thors SESSESHSSIH- Here, we relate 
the problem to similar bath models which have derived 
within the HAM formalism. 

The Hamiltonian Hs(t) implements the adiabatic 
quantum search, and will be written as 

H s (t) = /(t)(/-|* >< *o|)+<?(*)(I-|m >< m\) (27) 

where |\&o > corresponds to the initial state of the 
system, which we take as the equally-weighted superpo- 
sition l^o >= SiLi K > an d I is the identity opera- 
tor. The functions f(t) and g{t) will vary slowly during 
the running time to, and satisfy /(0) = 1, g(0) = , 
/(£g) = 0, g(tc) = 1- There are many possible choices 
of these functions, depending on the trade-off between 
time and energy cost one pursues [13, E3, US] ■ Here we 
choose these functions as obtained form imposing a local 
adiabatic condition [Vq . with f(t) = 1 — s(t) , g{t) = s(t) 
. In the large N limit 

s(t) = - H ^=tan(^= - arctanViV), (28) 

2 WN 

where e is a small number that controls the probability 
of success for the algorithm, which will run during a time 
t G = n^N/2e. 

The adiabatic quantum search evolution can be ef- 
fectively reduced to a two-level system, in the space 
spanned by the orthogonal vectors {\m >,\p >} , with 
\P > = y'jv-i J2i^ m N >• The minimum energy gap in 
the Hamiltonian lf27|) appears for eigenstates which are 
linear combinations of {\m >, \p >} . It occurs when 
s(t) ~ 1/2 and takes the value 1/y/N. The rest of eigen- 
vectors are degenerate, with eigenvalue f+g — 1, and are 
well separated from the previous two, specially around 
the avoided crossing point s ~ 1/2. It is reasonable to 
assume that one can restrict the system evolution to this 
effective two-level space. Arguments to favour this as- 
sumption are shown in [26]. 

We now introduce a model to describe the interaction 
with the environment. We will make an explicit com- 
parison with a numerical simulation of the Schrodinger 
equation. To this end, we make use of a simple model 
consisting on a band of Ni equally spaced levels. As we 
show, this model can describe relaxation to equilibrium 
and decoherence effects in a natural way, and may be 
regarded as a simplified version of the two-band model 
described in [ll[ . The Hamiltonian describing the envi- 
ronment is given by 

H * = Elf "I") ("I ( 29 ) 

n=l 1 



and the interaction between both systems by 

n 

V = ^]cr! h B l + h.c., (30) 

i=l 

with 

B l = \ i ^2 Cj(ni,n 2 )|ni)<n 2 |. (31) 

«2>>ll 

The indices n, ri\ and ri2 label the levels of the energy 
band, and o\ are Pauli matrices acting on qubit i . The 
global strength of the interaction with each one of the 
qubits is given by A^. The coupling constants Ci(ni,?i2) 
are independent Gaussian random variables. In order 
to make the model simpler, we will choose the same 
couplings for all qubits, which amounts to the replace- 
ment Y^i=i ^iCi(n%, n 2 ) — > n\c(ni,ri2) in Eq. (|30|) . The 
averages (denoted by <>) over the random constants 



c(ni,ri2) satisfy: 

(c(m,n 2 )) = 0, (32) 

(c(n 1 ,n 2 )c(n' 1 ,n' 2 )) = 0, (33) 

(c(n 1 ,n 2 )c*(n' 1 ,n' 2 )) = S^^S^^. (34) 



We will assume that an average over the possible real- 
izations of these coefficients is made when evaluating Eq. 

(ESI. 

Up to now, our model describes the coupling of the 
n qubits of the quantum computer to the environment. 
According to the above discussion, we will make the as- 
sumption that only the subspace spanned by the states 
{|m >, \p >} is relevant for the dynamics. Accordingly, 
we need to compute the matrix elements of Eq. lj30j) in 
this basis. A straightforward calculation gives, in the 
limit of large N = 2": 

V = a z C, (35) 
where a z acts on the system subspace, and 

n 

C=--J2(Bi + Bj) (36) 
i=i 

acts on system E. In the interaction picture, the above 
operator becomes 

C(t) = ~ c(n ll n 2 )e~ lt ^^\n 1 )(n 2 \ + h.c, 

Ji2>ni 

(37) 

with Lu(m, n 2 ) = j^(n 2 — n{) . 

The interaction Hamiltonian Eq. f35|) is of the general 
form l(20|) . with only one term appearing. Consequently, 
only one correlation function arises, which we represent 
by G(s, s') . Associated to this function, it exists one 
function T(t) defined as in (f24|) . One can also check that 
condition (|26|l is satisfied. 

As for the state pe , we make the simplest choice, 
by taking ps = jf- , where Ie is the identity operator 



in the environment space. This choice obviously satisfies 
[He, Pe) = . As a consequence, the correlation function 
G(s, s') only depends on the difference s — s' , and the 
function T(t) becomes independent of t . A redefinition 
of the variables s and s' 



r = lim — 

T->0 r 



gives 



ds / ds'G(s'). 
i] Jo 



(38) 



The correlation function can be calculated straightfor- 
wardly in the present model. We will consider the limit 
N% ^> 1 . In this limit, one obtains 



G(s)=N 1 ( 



71 A . &s. 2 

sin ) . 

26es 2 ; 



(39) 



The evaluation of the limit in Eq. (j38|) deserves some 
discussion. As will be shown below, we will be interested 
in time intervals which are much larger than 1/Ss . For 
such long-time variations, we can still consider values of 
t that are larger than l/5e, i.e. we assume that rSe ^> 1 
(see [13] for an extensive discussion). Within this con- 
text, Eq. (|38| finally gives 



n 2 X 2 TTN! 
8Se 



(40) 



The master equation ([23|) can be finally written, for 
our model, as 



dt 



Ps (t) = -i[H s (t),ps(t)}+r{<j z ps(t)<j z -ps{t)). (41) 



We have numerically solved Eq. (|41"Tl using the model 
presented above, for a system of n = 12 qubits. We 
choose a value of e = 0.1, for which the Grover time is 
to ~ 10 3 . The numerical values for the bath model 
are iVi = 2000 , Se = 0.5 . Therefore 1/Se = 2 < t G 
, in agreement with the approximations discussed previ- 
ously. We compare our numerical results to the solution 
of the Schrodinger equation of the total (S+E) system. 
The initial state is |*(0) >= |*s(0) > ®\V E (0) > , 

where |* s (0) >= j-\s > +\f^\p > and |* B (0) >= 



~Tkt S n ii \ n - > ' consistent with the above choice of pe 



According to the analysis based on the HAM equa- 
tions, one expects that Eq. (|4T]) will work when the con- 
ditions 



KffNi 1 

Cl ^^r~ * 2 

« 1 



(42) 



are met where the definition A e // = \n\ arises 

as a consequence of Eq. (|37|) . 

The results of our calculations are shown in Fig. 
1. We plot the probability of overlapping with the 
searched stated during the evolution of the system. Solid 
(blue) lines are obtained from the exact solution to the 
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Figure 1: (Color online) Numerical solution of the quantum 
search algorithm using the exact Schrodinger equation (blue, 
solid curves) or the derived master equation (J4TJ) (red, dashed 
curves). We plot the probability of overlapping with the 
searched state. Blue (solid) curves are obtained, from top 



to bottom, for A = , A = 10 and A 



5 x 10" 



For 



the two last cases, the corresponding approximated solution 
is also shown, also from top to bottom. Also for these two 
cases, green (dotted) curves show how decoherence manifests 
as time evolves, by plotting the magnitude defined in Eq. 

63). 



Schrodinger equation. The upper curve corresponds to 
the case of no coupling to the environment (A = 0) . 
Since the quantum computer runs during a 'Grover time', 
the probability approaches unity. The two lower curves 
have been obtained, from top to bottom, for A = 10~ 4 
and A = 5x 10~ 4 , giving Ci = 1.2, c 2 = 7.2 x 10 -4 and 
ci = 6,C2 = 1.8xl0~ 2 , respectively. The interaction with 
the environment translates into a worse performance of 
the adiabatic search, which is manifested as a lower prob- 
ability of success. This effect becomes stronger as the 
coupling to the bath increases. The degree of decoher- 
ence can be measured by several means. Here, as a figure 
of merit we calculate the magnitude 0|: 



C = j2Tr(pl) - 1, 



(43) 



which is also shown in the same figure for the same val- 
ues of A . Clearly, the decoherence increases with time. 
This effect is more pronounced for a larger coupling, giv- 
ing rise to an almost completely incoherent, and equally 
probable mixture, of the {|m >, \p >} states. 

One can also observe that the approximation obtained 
by solving the derived master equation (red, dashed 
curves) becomes more accurate for lower values of the 
coupling, in accordance to criteria Eq. lj42j) . Indeed, 
most of the difference observed for the smallest A are due 
to oscillations , corresponding to the fact that the com- 
plete numerical solution has been obtained for a partic- 
ular realization of the couplings in Eq. l[3"Tj). while these 
constants have been averaged out in obtaining (|4T1) . 
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VI. CONCLUSIONS 

In this paper, we have first discussed the so-called 
Hilbert space Average Method, as an alternative to de- 
scribe open quantum systems. We extended the method 
to the case of a system subject to a time-dependent 
Hamiltonian. We also made a connection of the evolution 
equations for this method with known master equations. 

We next discussed a simple model which can be use- 
ful for the study of a quantum computer performing an 
adiabatic quantum search, while in contact with an en- 
vironment. The ultimate purpose of such study is, of 
course, the understanding of the effects of decoherence 
on the performance of the computation. The model for 
the environment is simply a band of equally spaced lev- 
els with random coupling to the qubits of the quantum 
computer. In spite of its simplicity, we have shown that 
it incorporates decoherence effects in a clear way. The 
equations for the reduced system can be studied either 
under the form of HAM dynamics or master equations. 

One can also, for this model, perform an exact numeri- 



cal simulation of the total system (including the environ- 
ment). We have performed such a numerical study, and 
compared the results with the approximated dynamics of 
the system. As expected, increasing the strength of the 
coupling between the system and the environment im- 
plies a larger degree of decoherence, which translates into 
a lower probability of success for the quantum search. On 
the other hand, increasing the coupling also means that 
the master equation gives a poorer description of the ac- 
tual dynamics. The degree of approximation is controlled 
by the criteria derived for HAM equations within similar 
models for the environment. 
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